High energy neutron scattering from hydrogen using a direct geometry spectrometer 
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Deep inelastic neutron scattering experiments using indirect time-of-flight spectrometers have 
reported a smaller cross section for the hydrogen atom than expected from conventional scattering 
theory. Typically, at large momentum transfers, a deficit of 20-40% in the neutron scattering 
intensity has been measured and several theories have been developed to explain these results. 
We present a different approach to this problem by investigating the hydrogen cross section in 
polyethylene using the direct geometry time-of-flight spectrometer MARI with the incident energy 
fixed at a series of values ranging from E;=0.5 eV to 100 eV. These measurements span a much 
broader range in momentum than previous studies and with varying energy resolutions. We observe 
no momentum dependence to the cross section with an error of 4% and through a comparison with 
the scattering from metal foil standards measure the absolute bound cross section of the hydrogen 
atom to be a(H)= 80 ± 4 barns. These results are in agreement with conventional scattering theory 
but contrast with theories invoking quantum entanglement and neutron experiments supporting 
them. Our results also illustrate a unique use of direct geometry chopper instruments at high 
incident energies and demonstrate their capability for conducting high-energy spectroscopy. 

PACS numbers: 74.72.-h, 75.25.+Z, 75.40.Gb 



I. INTRODUCTION 

There has been great interest in the neutron scatter- 
ing cross section from hydrogen since the experiments of 
Chatzidimitriou-Dreismann et alX This initial work in- 
vestigated the scattering from hydrogen in mixtures of 
H2O and D2O and found an anomalously low value of 
the hydrogen cross section. The results are summarized 
in Fig. [1] where the ratio of the hydrogen and deuterium 
cross sections {(r{H) / a{D)) as a function of deuterium 
content is illustrated in panel b) and panel a) illustrates a 
representative time-of-flight spectra from which the cross 
sections were derived. The dashed line in Fig. Q] b) rep- 
resents the expected value of a(H) j a (D) =10. 7 based on 
conventional scattering theory. The data clearly deviate 
from expectations and were initially explained in terms 
of entanglement between hydrogen and deuterium atoms 
on very short (attosecond) timescales. 

Since this initial experiment, similar hydrogen cross 
section deficits have been observed in metal hydrides and 
hydrogen containing metallates (Refs. BIMJIals MMSl) , 
polystyrene and benzene (Ref. flfj ). Formvar (Ref. fill ) 
and HC1 (Ref. [l2r ). All of these results have been ob- 
tained with high-energy neutron scattering at large mo- 
mentum transfers using indirect geometry instruments 
(in particular the Vesuvio/eVs spectrometer at ISIS). 
The results, if correct, would represent a significant fail- 
ure of conventional models of scattering and several theo- 
ries have been developed to explain these results includ- 
ing the breakdown of the Born-Oppcnheimer approxi- 
mation (Ref ITsjQjilll) and quantum entanglement (Ref. 
0HM3I3) at short timescales. The latter theory has 
been refuted in Refs. [l9lr20lr2ll . where it has been shown 
to be inconsistent with the first-moment sum rule and 



also argued that good energy resolution and low tem- 
peratures are required to observe the effects suggested to 
result from entanglement. Such a sum rule violation docs 
not occur in the former theories where a redistribution 
of intensity conserves the first-moment while resulting in 
an apparent decrease in intensity of the hydrogen scat- 
tering. However, Ref. claims the former theories give 
too small an effect to be consistent with experiments. 

While supporting evidence has been claimed from elec- 
tron scattering (Ref. |H), Raman spectroscopy (Ref. I25T ). 
and neutron diffraction (Ref l24l) . several noteworthy at- 
tempts using neutron scattering have failed to reproduce 
a deficit in the cross section. Neutron interferometry has 
been conducted on mixtures of H2O and D2O and have 
found no deviations from expectations based on two flu- 
ids^ Several indirect geometry neutron scattering mea- 
surements utilizing the foil filter technique, an experi- 
mental setup similar to the initial study presented in 
Ref. [l|, have been conducted and report no evidence for a 
missing hydrogen cross sectioni 27 i 28 ' 29 These studies have 
been reconciled with the initial data based on the large 
difference in the energy resolution and hence have ques- 
tioned whether da/dfl or d 2 /dEd£l was being measured 
and compared i 30 i 31 ' 32 A detailed revisiting of the initial 
experimental setup on polyethyelene did find evidence for 
an anomalous lack of intensity consistent with previous 
results, however several experimental concerns were also 
notedi^ 3 . 

Given the large discrepancy in experimental results 
and the theoretical interest in this effect, we have pur- 
sued a different approach to the investigation of the hy- 
drogen cross section as a function of momentum transfer. 
While all previous experimental studies have been con- 
ducted with indirect geometry spectrometers (where the 
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final neutron energy is fixed at one energy), we present 
a direct geometry study of the cross section of hydrogen 
in polyethylene using the MARI chopper spectrometer. 
Because the incident beam energy can be fixed at a par- 
ticular value, the hydrogen cross section can be inves- 
tigated over a range of momentum transfer with vary- 
ing energy resolution and this allows a direct comparison 
with the experiments and theories described above. This 
experiment also allows overlapping data to be obtained 
in both momentum and energy transfer thereby checking 
the consistency of the experimental results. This feature 
of direct geometry spectrometers is not possible on indi- 
rect machines where a filter edge or analyzer forces one 
particular fixed final energy. 



Given the previous results on all hydrogen containing 
materials, there are two key results which need to be 
investigated and provide the basis for the experiment re- 
ported here. Firstly, the momentum dependence of the 
scattering line shape needs to be measured over a broad 
range extending to over Q ~ 100A -1 . Secondly, the ab- 
solute value (in units of barns) is required for these mo- 
mentum transfers and comparison with the theoretical 
value made. 



We will show that within experimental error, there is 
no change in the hydrogen cross section with momentum 
transfer and that the absolute value of the cross section is 
consistent with conventional scattering theory. These re- 
sults show that conventional scattering theory based on 
the Born-Oppenheimer approximation and the impulse 
approximation is valid and provides a good description 
of the data. These results also demonstrate the use of di- 
rect chopper instruments for conducting spectroscopy at 
high energies and we suggest possible future uses of these 
instruments for understanding problems in magnetic and 
strongly correlated electronic systems. 



The paper is divided into five sections and two appen- 
dices. The first section discusses the experimental setup 
including the direct geometry MARI spectrometer and 
also how we have measured and subtracted the back- 
ground. The second section reviews the impulse approx- 
imation and the data analysis which provides the basis 
for understanding the experimental results. We also re- 
view the kinematics of a direct geometry spectrometer. 
The third section describes the calibration of the spec- 
trometer using a series of "standard" metal foils such 
that the experiments could be put on an absolute scale. 
The fourth section describes the hydrogen cross section 
and outlines the momentum dependence and the abso- 
lute value obtained from the comparison with the results 
from the metal foils. We end the paper with a section 
of discussion, conclusions and future plans. The appen- 
dices include discussions on the energy calibration and 
resolution. 
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FIG. 1: The initial results illustrating an anomalous deficit of 
intensity in H2O and D2O. The time-of-flight signal is plotted 
in panel a) which illustrates the distinct separation observed 
between the H and D recoil peaks. The ratio of the intensities 
as a function of D2O concentration is shown in panel b. The 
dashed line at 10.7 is the expected value based on conventional 
scattering theory. The results are taken from Ref. QJ. 



II. EXPERIMENT 

A. The MARI direct geometry spectrometer 

MARI is a direct geometry time-of-flight spectrometer 
as shown in Fig. [2] a). The spectrometer consists of a 
Fermi chopper 10.02 m from the neutron moderator, a 
sample position 1.77 m from the chopper and a detector 
bank of He 3 tubes located 4.02 m underneath the sam- 
ple. The incident energy is determined by the phasing 
of the Fermi chopper with respect to the neutron pulse 
produced in the moderator while the scattered energy 
is determined by the time-of-flight of the scattered neu- 
trons. 

For energy transfers below ~ 1 eV, both the nimonic 
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FIG. 2: A schematic of the direct geometry (a) MARI and 
indirect geometry (6) VESUVIO spectrometer, a) The lo- 
cations of the beam monitors (labeled Mon 1 etc.), nimonic 
and Fermi choppers, and detectors are marked. The incident 
beam energy is fixed by the Fermi chopper and scattering neu- 
trons are detected by the He 3 detectors, b) No Fermi chopper 
exists on VESUVIO and an incident beam with a spectrum 
of energies is incident on the sample. The final energy is fixed 
by using an Au foil resonance. 



and Fermi chopper are utlized. The nimonic chopper re- 
moves high-energy neutrons (typically above ~ 3-5 eV) 
and other high-energy radiation that is emitted from the 
moderator at very early times after the arrival of the pro- 
ton pulse on the target. For our high energy neutron ex- 
periments with incident energies above 1 eV, the nimonic 
chopper was removed so that the high-energy spectrum 
could be accessed. This inevitably resulted in a signif- 
icant background at short times on the detectors. The 
subtraction of this background and spectrometer calibra- 
tion are discussed later. 

The Fermi chopper consisted of blades of neutron ab- 
sorbing boron nitride. The blades were curved such that 
for incident energies of 500 meV, the neutrons would have 
an unobstructed path when the chopper is spinning at 
600 Hz. Various other blade configurations are available 
for lower incident energies. The incident energy was fixed 
by adjusting the phase between the opening of the Fermi 
chopper and the time the proton pulse reached the target. 
The incident energy was further calibrated for the effects 
of drift in the detectors and possible timing errors by 
recording the time the neutron pulse reached the Mon 2 
and Mon 3 glass bead detectors. This calibration is out- 
lined in detail in Appendix A to this paper. Because the 
distances are all known, a time-of-flight signal in both 
these monitors provides an unambiguous measurement 
of the incident neutron energy. Mon 1 was used to nor- 
malise all data sets to an equal incident beam flux by 
dividing by the integrated counts over all times-of- flight. 

The detectors were cylindrical tubes of He 3 at 10 atm 
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FIG. 3: The scattered intensity off a vanadium foil with a 
white beam of neutrons incident on the sample as measured 
on the detectors. The dashed line is with the nimonic chopper 
in and the solid line is with the nimonic removed. The data 
is integrated for 29 = [5°, 130°]. 



of pressure. The detector efficiency as a function of en- 
ergy was calculated based on the known cross-section 
of He 3 assuming a cylindrical geometry^ Counts were 
recorded on the detector with a constant time intervals 
of (5i=0.25 /is up to 1000 /is from when the proton pulse 
is incident on the target. To convert counts measured on 
the detectors to absolute units (barns) we used a series of 
metal foils with known tabulated cross sections. We de- 
scribe the method of obtaining the calibration constant 
later in this paper. 

We conducted two different experiments. The first 
measured the scattering of polyethylene wrapped around 
an aluminum annulus. Given the need to obtain absolute 
cross sections for the hydrogen atom and difficulties in re- 
moving the aluminum intensities, we performed a second 
experiment using samples of either low-density (LDPE) 
or high-density (HDPE) polyethylene (CH 2 ) films. The 
films and metal foils were free standing and mounted at 
42° to the vertical axis. This setup had the advantage 
that only the sample was in the beam for any given mea- 
surement and this simplified the analysis considerably 
without the need to consider the presence of aluminum 
in the scattered signal. Both experiments gave consistent 
results and we focus on the second experiment based on 
free standing films in this paper. The physical properties 
(including thickness d, density p, and molar mass m) of 
the polyethylene samples are listed in Table [IJ 
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TABLE I: Polyethylene Samples 



X 


d (mm) 


P (g/cm 3 ) 


m (g/mol) 


HDPE 


0.044 


0.95 


14.03 


LDPE1 


0.20 


0.92 


14.03 


LDPE2 


0.40 


0.92 


14.03 



Lead 



It is interesting to compare the MARI direct geometry 
spectrometer described here with VESUVIO which was 
used for the earlier experiments on hydrogen. VESUVIO 
is an inverse time-of-flight machine which has a pulsed 
white incident beam of neutrons that are scattered from 
the sample and the scattered energy is determined by the 
transmission through a thin film for which there is a nu- 
clear resonance such as gold. The energy of the incident 
neutrons is known from the time-of-flight of the neutrons 
and energy of the scattered neutron is determined by 
the nuclear resonance. This spectrometer has the disad- 
vantage that it involves the subtraction of counts with 
the gold foil 'in' and the foil 'out'. It also depends cru- 
cially on the energy dependence of the incident neutron 
flux. Since there are very few metal foils that have suit- 
able resonance energies and widths for use in this type 
of spectrometer (outlined in Ref. [H]) it is not possible 
to vary appreciably the scattered neutron energy and so 
test the instrumental results by measuring a particular 
wavevector and energy transfer with different neutron 
energies. In contrast, a direct time-of-flight instrument 
such as MARI has a monochromatic incident beam and 
the detectors measure neutrons with a range of different 
energies. The resultant flux is then crucially dependent 
on the knowledge of the efficiency of the detectors as a 
function of the neutron energy. The incident energy for 
this type of instrument can, however, be controlled and 
varied and this enables experiments to be performed to 
ensure that measurements at different incident energies 
give consistent results. 



B. Background characterisation and subtraction 

As described in the previous section, because of the re- 
quirement of large incident energies the nimonic chopper 
was removed for these studies. This had the effect that 
(while allowing access to high-energy neutrons) it also al- 
lowed high-energy charged particles and 7-rays from the 
initial interaction of the proton pulse with the tungsten 
target to be incident on the sample. This contributed 
to a significant increase in the background especially at 
short times of flight. This is shown in Fig. [3] which il- 
lustrates the counts measured on the detectors integrated 
over 29 = [5°, 130°] with the nimonic chopper in and out. 
The Fermi chopper was removed for these measurements 
to ensure a white beam spectrum incident on the sample 
position. The sample was a vanadium foil of thickness 
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FIG. 4: Time of flight spectra for a lead foil sample for 
Ei=2.5 eV, 40 eV, and 100 eV. The data is integrated over 
2(9=[30°,70°]. to represents the time the pulse of neutrons 
strikes the sample and is graphically represented by the ver- 
tical dahsed lines. 



0.075 mm. 

Fig. [3] illustrates that at small times, below ~ 1000 fis, 
the nimonic chopper has a significant effect on the spec- 
trum and counts on the detector. A particuarly strong 
feature is a significant tail which seems to diverge as t ap- 
proaches t = 0. It is essential to understand this strong 
background scattering so that it can be quantitatively 
subtracted to obtain the correct intensity and line shape 
of the scattering. We found that the background is to 
a good approximation independent of the incident neu- 
tron energy chosen. This is demonstrated in Fig. Uwhich 
shows the scattering for several different incident energies 
from a lead foil. The dotted vertical line shows the time 
when the neutron pulse reached the sample (to) and so for 
earlier times, there can be no effect of the incident energy. 
Hence, then by performing scans at different energies the 
background spectra can be constructed as shown for lead 
foil in Fig. [5) The result is then independent of the inci- 
dent neutron energy although it depends on the sample so 



5 




50 100 150 200 250 300 



tfrlS) 

FIG. 5: The time-of-flight spectra for lead foil plotted for 
Ei=20, 40, and 100 eV. The background derived from the 
method described in the text is plotted as the solid and thick 
black line. 

a similar procedure was followed for each sample studied. 
Representative plots of the background subtracted time- 
of-flight spectra for lead are plotted in Fig. [HI The plots 
show that the background subtraction described above 
does accurately remove the background intensity as the 
data subtracts to zero for both shorter and longer times 
than the lead recoil peak. We describe in the next sec- 
tion how this time-of-flight spectra was converted into an 
energy spectrum. 
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FIG. 6: The time-of-flight spectra with a lead foil placed at 
the sample position and with the background subtracted by 
the method described in the text. The intensities were ob- 
tained by integrating 29— [30,70]. The background is zero 
within the errors at times both shorter and larger than the 
peaks in the spectra. 



III. KINEMATICS AND THE IMPULSE 
APPROXIMATION 

In this section we review the kinematic equations based 
on energy and momentum conservation and line shape 
predicted by the impulse approximation. Having dis- 
cussed how we obtained a background subtracted data 
set in time in the previous section, we will also outline 
the conversion to energy spectra which can be compared 
directly with theory. 

A. Kinematics and equations for detected neutron 
intensity 

Based on conservation of energy and momentum, the 
peak positions and the neutron cross sections can be cal- 
culated. If we consider a beam of incident neutrons with 
a wavevector fco which scatters off a sample of hydro- 
gen atoms resulting in a scattered wavevector of k\ , the 
energy transfered is related to the momentum of the re- 
coiled hydrogen atom by E — h 2 Q 2 /2M p , where M p is 



the mass of the proton. Making the approximation that 
the neutron has the same mass as the hydrogen atom, en- 
ergy conservation results in Q 2 = k 2 — k 2 . Taking conser- 
vation of momentum gives Q 2 = k 2 + k\ — 2k$k\ cos(26*), 
with 26 defined as the angle between the incident and 
final neutron wavevector or the scattering angle. Sub- 
tracting the energy and momentum conservation equa- 
tions results in 

k x = fc cos(26>), (1) 

while adding gives Q — fcosin(2#). This latter equa- 
tion illustrates that the maximum momentum and energy 
transfer occur when 20=90°. 

The measured intensity of the scattered neutrons as a 
function of 29 needs to be corrected for several factors. 
The scattered cross section per unit incident flux into 
solid angle d£l and an energy interval dE (d 2 a/d£ldE) 
for a monotonic system is proportional to ki/koS(Q, E) 
where S(Q, E) is the van Hove scattering function. The 
derivation of S(Q, E) assumes an incident flux of neu- 
trons falling on the sample. The equations and correc- 
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tion factors for an indirect geometry machine (such as 
VESUVIO) have been described in RefsQl andlM Here 
we outline the appropriate equations and corrections for 
a direct geometry spectrometer. 

The flux on the sample depends on the moderator and 
the properties of the Fermi chopper. The flux from the 
moderator is proportional to C(Eo)dEo where C(Eq) is 
approximately proportional to 1/Eq. While it is impor- 
tant to know the flux as a function of energy on a indi- 
rect geometry spectrometer where only the final energy 
is fixed, the details of the incident beam spectrum are 
not crucial to the use of a direct geometry spectrometer 
where the incident flux and energy are fixed through the 
phasing and frequency of the Fermi chopper. 

Since the chopper spins at a constant speed the flux 
coming through the chopper is proportional to the open- 
ing time which is almost independent of the phasing 
and hence of the energy. The chopper then lets neu- 
trons through for a constant time interval and so the 
flux through the chopper is C(Eo)dEo/dtAt where At 
is the chopper opening time. Since E = l/2mv 2 = 
l/2m(£/t) 2 , dEo/dt — m£ 2 /t 3 and hence is proportional 

3/2 

to E . These equations combined give the following for 
the cross section, 
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where A is a constant of proportionality. 

The detectors in a direct geometry spectrometer record 
intensity as a function of time and hence measure 
d 2 a/dftdt rather than d 2 a/dildE. The conversion fac- 
tor between the two is dE\jdt and has been calculated 

3 /2 

above to be proportional E ± ' . These relations can be 
put together to give the following expression for the mea- 
sured cross section 
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Setting Eqn. [2] equal to Eqn. [3] gives the following ex- 
pression. 



FIG. 7: The intensity as a function of energy transfer for a 
lead foil with Ei=20, 40, and 100 eV. The data is the same 
as presented in Fig. [6] though converted to an energy spec- 
trum and corrected for the detector efficiency. The data was 
obtained by integrating over 29— [30,70]. 



S{Q,E) = BE 2 



dVldt 



= B't 4 
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(4) 



where B and B' are constants. The expression assumes 
a direct geometry spectrometer with a fixed incident en- 
ergy and a fixed phasing of the incident beam Fermi 
chopper. The parameter t is the time taken from when 
the neutron pulse hits the sample and when a signal is 
recorded on the detectors, ft should be noted that the pa- 
rameter t is the time from when the proton pulse strikes 
the target and these two different definitions of time will 
be maintained throughout this paper. Eqn. 0] illustrates 
that S(Q, E) is related to the measured intensity (cor- 
rected for detector efficiency) after a correction of t 4 . An 



example of the background subtracted data converted to 
an energy spectra is illustrated in Fig. [7]for the data pre- 
sented previously in Fig. [H] for a lead foil. We note that 
the energy dependent detector efficiency was corrected 
for when the data was converted from time to energy. 

The width and the integrated intensity depends on the 
trajectory in the E — Q plane of the measurement. So 
far in the discussion of the width we have assumed that 
we are using a constant wave vector trajectory and are 
varying only the energy. Since the data is collected by 
detectors at certain scattering angles it can be advanta- 
geous to use a constant scattering angle trajectory when 
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problems arising from changing detectors across a scan 
do not arise. The disadvantage usually is that it is then 
necessary to use a Jacobian as first introduced by Waller 
and Froman^i For a conventional time-of-flight machine 
with hydrogen as the sample the Jacobian at the peak of 
the scattering is 



J(2fl) = l + g(l-gco B (2(9)). 



(5) 



where m is the mass of the neutron and M is the mass 
of the sample atom. For the case of hydrogen, m=M 
and using this result in conjunction with the kinemati- 
cal scattering expressions (Eqn. [T]) we obtain that the 
Jacobian is identically unity. Consequently to a good 
approximation the intensity and width can be found di- 
rectly by integrating the spectra at each scattering angle. 
The width is then the same as that obtained for a con- 
stant Q trajectory and the intensity has the same value 
as for a constant Q trajectory. 



B. Impulse approximation and line shape 

To understand the observed cross sections and scatter- 
ing profiles, we use the impulse approximation which is 
the high energy neutron scattering from a single atom av- 
eraged over all the different atoms in the sample. The in- 
terference terms in the scattering amplitude are assumed 
to completely cancel. The scattering of neutrons with an 
incident energy of Eq and wave vector fco to a state with 
wave vector k± and energy E\ is given by 



1(E) = Aexp - 



(E-CQ 2 f 
(WQf 



(6) 



Where E is the energy transfer Eq — E\ and Q is the 
wave- vector transfer Q = ko — hi. The scattering func- 
tion is normalised if A = (^/^WQ)^ 1 and the constant 
C = fi 2 /2M where we shall take M as the mass of the 
hydrogen atom. The width of the scattering arises from 
the motion of the hydrogen atom immediately before the 
neutron collides with it and W is proportional to the 
mean average velocity of the hydrogen atoms. Eqn. [6] 
demonstrates that the width is proportional to momen- 
tum transfer Q and hence the width of the recoil line is 
predicted to increase with increasing scattering angle and 
hence momentum transfer. 

It is useful to introduce dimensionless coordinates to 
describe the scattering by dividing the energies by the 
incident neutron energy Eq and the wave vectors by the 
magnitude of the incident neutron wave vector fco- The 
expression for the scattering then becomes 



1(e) =Bexp - 



( e~q 2 
(wq) 2 



(7) 
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FIG. 8: The measured integrated intensity is plotted as a 
function of the calculated intensity per unit area for a sc- 
ries of calibration foils. The thicknesses and cross sections 
are discussed in the text. The intensities were obtained by 
integrating 20=[3O,7O]. 



Where the reduced incident energy e = E/Eq and q — 
Q/ko. The normalisation constant is B = A/ko and 
the reduced width is w = W/(Cko). This expression 
shows that the peak energies for the scattering occur 
when e = q 2 and that this is independent of the incident 
neutron energy showing that a plot of the peak energies 
should be identical for all incident energies when plotted 
in terms of the reduced coordinates, e and q. Rearrang- 
ing the neutron scattering expressions for conservation of 
energy and momentum and using the expression for the 
peak of the scattering the wave vector transfer becomes 
q = sin(20) where 29 is the scattering angle. Again, this 
equation illustrates that the momentum transfer q is al- 
ways less than 90°. 

The half-width of the scattering in these reduced co- 
ordinates is given in terms of the reduced energy as 
T(q) = qwVhi2. This depends on the angle of scatter- 
ing q and is also proportional to the incident wave vector 
fco. Hence as the incident energy increases w decreases 
and the width in a reduced plot also decreases. Our con- 
clusion is that when using these reduced coordinates the 
position of the peak should be independent of the inci- 
dent energy but that the peak in the scattering should 
steadily get narrower as the incident energy increases. 
Furthermore, if polyethylene can be approximated as a 
harmonic solid the line shape of the scattering is Gaus- 
sian and therefore we use this linsehape to investigate the 
recoil scattering in detail. 
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FIG. 9: The scattering from a 0.4 mm thick film of LDPE as 
a function of time of flight. Panels a) — c) represent constant 
29 cuts at 5°, 25°, and 50° respectively. Panel d) is a false 
color plot of the intensity as a function of 29 and E. 



IV. INTENSITY CALIBRATION AND 
ABSOLUTE UNITS 

To convert the measured intensities to an absolute 
cross section with units of barns, a series of standard 
materials with known cross-sections were compared. We 
have chosen to use foils of copper, lead, nickel, niobium, 
and vanadium. The thicknesses of the metal foils are 
listed in Table ITTl and were chosen to scatter less than 10 
% of the incident beam so that corrections due to multi- 
ple scattering would not need to be considered. The use 
of heavy atom foils with masses much larger than the 
hydrogen atom means that the recoil energy was small 
in comparison to both the recoil energy of the hydrogen 
atom and the experimental resolution. 

The cross sections, a were taken from the tables com- 
piled by Sears (Ref. [36|) for neutron energies of 25.3 meV. 
The metal foils chosen have cross sections which vary lit- 
tle over the energy range of interest and do not display 
strong resonances below ~ 100 eV. 



TABLE II: Foil Standards 

X d (mm) p (g/cm 3 ) a (barns) m (g/mol) 



V 


0.085 


6.11 


5.1 


50.94 


Ni 


0.11 


8.908 


18.5 


58.69 


Nb 


0.51 


8.57 


6.3 


92.91 


Pb 


0.86 


11.34 


11.1 


207.2 


Cu 


0.73 


8.94 


8.0 


63.55 



From the values in Table [TT1 we calculated an expected 
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FIG. 10: The scattering from a 0.4 mm thick film of LDPE as 
a function of energy transfer. Panels a) — c) represent constant 
29 cuts at 5°, 25°, and 50° respectively. Panel d) is a false 
color plot of the intensity as a function of 29 and E. 



cross section per unit area for each foil (in units of 
barns x mol/cm 2 ) equal to I ca i — dpa/m where d is 
the thickness, p is the density, and m is the molar mass. 
The measured intensity (integrated over 2f?=[30°,70°]) 
was then plotted against the calculated cross section to 
obtain a calibration constant otfoils from the slope. The 
calibration curves for Ei=20, 40, and 100 eV are illus- 
trated in Fig. [5] The measured intensity scales linearly 
with the calculated intensity and therefore confirms that 
multiple scattering effects do not need to be considered. 
The systematic increase in slope with incident energy rep- 
resents the fact that the flux of neutrons increases with 
incident energy. 

The curves provide an accurate means of putting the 
measured integrated intensity on an absolute scale. The 
increased scatter in the data points for higher incident 
energies is possibly due to uncertainty in the energy de- 
pendence of the cross section^ While the energy depen- 
dence of the cross sections have been tabulated, we have 
found substantial variations in the tabulated values (up 
to ~ 10%) and therefore have chosen to take the val- 
ues at 25.3 meV which have been measured and carefully 
tabulated. These values should therefore provide a good 
approximation of the bound cross section. Based on the 
slopes {a f u s ) obtained from Fig. [5J the integrated in- 
tensity can be converted to an absolute cross section in 
units of barns. We consider that a large component of the 
error is due to uncertainty in the energy dependence of 
the cross-sections in the range 1-100 eV. Improved mea- 
surements are needed to obtain more accurate results. 



9 



Ei=20 eV 



LDPE (0.4 mm) 

Ei=40 eV 



Ei=100 eV 




10 20 -20 



20 

E(eV) 



25 25 50 75 



100 



80 



> 

■Si- 



60 



40 



20 



LDPE (0.4 mm) 

□ 100 eV 
40 eV 
o 20 eV 



2 3 4 

Q 2 (1 4 Angstroms" 2 ) 



FIG. 11: The measured intensity plotted as a function of 
energy transfer with a sample of LDPE (0.4 mm). Panels 
a)— c) are for Ei=20 eV,d) — /) are for Ei=40 eV, g)—i) are for 
Ei = 100 eV. The data is plotted for representative scattering 
angles of 5°, 25°, and 50°. 



FIG. 12: The position of the hydrogen recoil line is plotted 
as a function of Q 2 for E;=20, 40, and 100 eV. The dashed 
line is E — h 2 Q 2 /2M P . The errorbars represent the full width 
at half maximum intensity of the hydrogen recoil peak. The 
fitted error bars are approximately the size of the data points. 



V. HYDROGEN SCATTERING 

In the previous sections we outlined the background 
subtraction, conversion of the time spectra to energy, 
and the calibration of the spectrometer using a series of 
metal foils. In this section we now discuss the results for 
polyethylene in both units of time-of-fiight and energy 
transfer. We then discuss the momentum dependence of 
the hydrogen recoil line line shape including the recoil 
peak position, width and integrated intensity. 



A. Hydrogen recoil in polyethylene 

Representative time-of-fiight spectra for the LDPE 
sample (d=0A mm) are illustrated in Fig. [9l The data 
in panels a — c) were obtained with an incident energy 
of Ei=40 eV at three different scattering angles (5°, 25°, 
and 50° respectively) sampled at an equal time intervals 
of 0.25 fis. Panel d) illustrates a false contour plot of the 
intensity as a function of both time-of-fiight and scat- 
tering angle. The energy converted data is displayed in 
Fig. [TU] and is based on the same data set as in Fig. 
[9] using the procedure outlined in the previous sections. 
Large energy transfers corresponds to longer time and 
since there are more time bins per unit energy transfer, 
the density of points increases with energy transfer. 

At the lowest scattering angle (28—5° in panel a) of 



Fig. [9] and Fig. IT0|) . the momentum transfer is suffi- 
ciently small that the recoil lines of hydrogen and car- 
bon overlap and appear at the elastic (E—0 in Fig. [TU]) 
position (appearing at a time of ~ 180 /is in Fig. [5]). 
Since the incident energy is fixed on MARI, increasing 
energy transfer correponds to larger times on the plots 
illustrated in contrast to indirect geometry spectrome- 
ters where energy transfer grows with decreasing time- 
of-fiight. 

As the scattering angle increases along with the mo- 
mentum transfer, we observe a gradual shift of the recoil 
line (20=25°) and then finally a separation of the car- 
bon and hydrogen recoil lines clearly observed at 20=50° 
with the carbon near E=0 and the hydrogen recoil peak 
appearing at finite energy transfer. Because the energy 
dependence of the recoil line for an atom of mass m is 
E = ti 2 Q 2 /2m, the carbon atom will recoil with a smaller 
energy than the hydrogen atom for a given momentum 
transfer. The recoil peak of carbon will therefore appear 
at smaller energy transfers, or at shorter times given the 
geometry of our time-of-fiight instrument. 

A series of different representative constant 20 cuts for 
the three different incident energies studied is illustrated 
in Fig. [TT] for scattering angles of 28 — 5°, 25°, and 
50°. All three different incident energies display the same 
qualitative features with a strong hydrogen recoil peak 
which shifts in energy as a function of 28 and a very 
weak carbon recoil peak located near the elastic position 
at all scattering angles. 
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FIG. 13: The full-width as a function of momentum transfer 
is plotted for Ei=20, 40, and 100 eV. The data is from a 0.4 
mm thick film of LDPE. 



B. Variation of Scattering with Q 
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The momentum dependence of the hydrogen recoil line 
was obtained by fitting the constant 29 scans to the sum 
of two Gaussians to represent the recoil lines from hydro- 
gen and carbon with a total of 6 parameters (2 ampli- 
tudes, linewidths, and positions). To restrain the fit and 
obtain an accurate measure of the hydrogen recoil posi- 
tion at small scattering angles, where both the carbon 
and hydrogen recoil lines overlap, we fixed the amplitude 
and linewidth of the carbon line to be that measured at 
scattering angles greater than 100°. Because the neutron 
and proton masses are approximately equal, there is no 
hydrogen recoil above a scattering angle of 90°. With 
these constraints, we were able to obtain a good descrip- 
tion of the line shape and scattering angle dependence 
for the three different incident energies of the neutrons. 

The peak position of the hydrogen recoil line is dis- 
played in Fig. [T2] with the errorbars representing the 
measured width for the different incident neutron ener- 
gies used. According to the impulse approximation the 
peak position should scale with momentum transfer Q as 
E = h 2 Q 2 /2M p and is represented by the dashed line in 
Fig. 1121 The agreement between experiment and theory 
is excellent over the entire range of Q and for the dif- 
ferent incident energies. We observe no deviation which 
would be expected if there was a breakdown of the Born 
Oppenhemier approximation predicted by several theo- 
ries^ 

The full-width of the scattering (2r) is shown in Fig. 
1131 The results in contrast to expectations based on the 
impulse approximation which predict an increase in the 



FIG. 14: The integrated intensity of the hydrogen recoil line 
as a function of 26 for a) E 4 =20eV and b) Ei = 100 eV. 



linewidth with increasing momentum transfer shows a de- 
crease. This figure makes clear that the width is largely 
dominated by the experimental resolution which, par- 
ticularly for incident energies of 100 eV, becomes much 
smaller at larger Q. As discussed in Appendix B, this 
is because it is dominated by the width of the incident 
neutron pulse from the target and moderator. With 20 
eV incident neutrons this is less of a problem and there 
is a slight increase in the full- width with increasing Q. 
At low angles the full width of the resolution is approxi- 
mately 15 %, 35 %, and 65 % for incident energies of 20 
eV, 40 eV, and 100 eV respectively and at large angles is 
25 %, 17 %, and 14 % at high angles. 

The integrated intensity as a function of scattering an- 
gle is plotted in Fig. [14] for the LDPE (d=0A mm) for 
Ej=20 (panel a)) and 100 eV (panel b)). The data sets 
represent the two extreme ranges of momentum trans- 
fer investigated and illustrate no measurable momentum 
dependence to the integrated intensity of the hydrogen 
recoil line. A summary of all of the data taken at each 
incident energy as a function of momentum transfer (Q) 
is plotted in Fig. EH The data in both Figs. [Hand US] 
have been obtained by dividing by the intensity for each 
incident neutron energy and sample extrapolated to its 
value at Q=0. The density of data points at lower mo- 
mentum transfers in Fig. E5]is larger owing to the signif- 
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FIG. 15: a) The intensity normalised as Q — > is plotted 
as a function of momentum transfer for Ei=20, 40, and 100 
eV. The mean of the data is 1.00 with a standard deviation 
of 0.04. b) A histogram of the of the observed intensities is 
plotted illustrating the fluctuations of the data. 



icant overlap of data points from each incident energy at 
low scattering angles. For the largest momentum trans- 
fers, only the Ei=100 eV data set contributes. Based on 
the data in panel a) we find no momentum dependence 
with the dashed line equal to 



HQ) 

I(Q = 0) 



= 1.00 ±0.04. 



(8) 



The error is derived from the standard deviation of the 
data points in Fig. [15] a). Panel 6) is a histrogram of 
the points in panel a) and visually illustrates the statis- 
tical spread of the data. Based on the combined data 
in Fig. 1151 we conclude there is no momentum depen- 
dence to hydrogen cross section. This is in contrast with 
previous results taken on indirect geometry machines (as 
illustrated in Fig. [1} and most notably those reported in 
Ref. which spanned a momentum transfer of up to ~ 
100 A" 1 . 



C. Absolute cross section 

Having shown that the intensity of the hydrogen recoil 
scattering is constant with momentum transfers ranging 
up ~ 250 A -1 , we now discuss the absolute value for the 
hydrogen cross section. We have obtained 9 independent 
measurements (three different densities at three different 
incident energies) of the recoil cross section in polyethy- 
lene and each measurement can be put on an absolute 



FIG. 16: The measured intensity as a function of calculated 
intensity normalised by the number of hydorgen atoms. 



scale through the calibration curves obtained from the 
measurement of metal foils described above. 

Figure [TH] summarizes the results for the three differ- 
ent polyethyelene samples measured at the three different 
incident energies. The calculated intensity is normalised 
by the hydrogen cross section in each formula unity of 
polyethylene (CH2) and therefore represents a measure 
of the number of hydrogen atoms per unit surface area. 
The measured intensities for each density and thickness 
vary linearly with the density and therefore illustrate that 
multiple scattering processes are negligible. 

The momentum dependence was found to be indepen- 
dent of the scattering angle for each density of polyethye- 
lene at each incident energy. The average intensity from 
each data set gives a different measurement of the re- 
coil cross section and there are 3 polyethylene samples 
and 3 different incident neutron energies giving 9 differ- 
ent measurements of the average intensity. The inten- 
sity of each measurement (I m es(H 2 ) m units of 'Counts') 
was plotted against the value af a s x I ca ujj 2 \/(r{H2) (in 
units of Counts/barns) represented in Fig. [TB] and the 
result is showin in Fig. [T7J The straight line has a slope 
which is the hydrogen cross section in one formula unit 
of polyethyelene denoted <j{H2). 

The results are displayed in Fig. [17] and illustrate that 
all of the data do fall on a straight line illustrating the 
consistency of each data set obtained with different inci- 
dent energies and differing energy resolutions. The slope 
of the line divided by 2 is the bound cross section of a 
single hydrogen atom, 



a(H) = 80 ± 4 barns, 



(9) 
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FIG. 17: A plot of the measured intensity as a function of 
the calculated intensity as described in the text. All three 
different incident energies and the three thickness of HDPE 
and LDPE are plotted and denoted by different symbols. The 
slope of the curve is the absolute cross section of a{H2)- 

in perfect a gree ment with value of 82.02 ± 0.06 quoted by 
Sears (Ref. [3fJ) for the bound cross section at an energy 
of 25.3 meV. The error bars in Fig. [TTland the resulting 
error bar on the final absolute cross section is the result 
of uncertainties in the foil thickness, foil cross section at 
high energies, and also the thickness and hence density 
of polyethyelene films. Nevertheless, our data is in good 
agreement with the expected cross section based on the 
impulse approximation. 

In comparison, previous results based on indirect spec- 
trometers have given values ranging from cr{H) =48-64 
barns at large momentum transfers of about 100 A -1 . 
The experiment discussed here clearly rules this value 
out over a much broader range in momentum transfer. 

VI. CONCLUSIONS AND DISCUSSION 

The experiment described above has shown conclu- 
sively three key results. Firstly, a direct geometry time- 
of-flight machine can be used with incident energies up to 
at least 100 eV and the background quantitatively sub- 
tracted. Secondly, the hydrogen recoil cross section is 
constant over the entire momentum transfer range of the 
experiment (Q ~ 250A -1 ). Thirdly, the absolute value 
of the cross section is in perfect agreement with previ- 
ous measurements of the bound cross section as expected 



from conventional scattering theory. All of these results 
are in contrast to the results discussed in the introduction 
on a variety of hydrogen containing systems. 

While we have only conducted experiments on 
polyethylene, the results are much more general. The im- 
pulse approximation assumes that the scattering sites are 
independent and therefore the results are indeed general 
to all hydrogen containing systems. Hence, our results 
stand in contrast to the variety of hydrogen containing 
systems where a deficit in intensity was observed. 

It is also of interest to discuss and compare the en- 
ergy widths of our data and those obtained using Vesu- 
vio when the cross section was anomalous. In Fig. I13l wc 
have presented all of our data in terms of energy transfer, 
while previous results were presented in terms of time. 
The expression relating energy resolution to the timing 
resolution can be obtained by considering a fixed scatter- 
ing angle 29 and a spectrometer configuration in which 
E = C/t 2 . Taking derivatives of this expression we ob- 
tain AE/E ~ 2At/t. Using the results of Fig. [TBI the 
energy resolution on MARI at large scattering angles is 
AE/E ~ 20%. The experimental data shown in Fig. 
UJ gives an energy resolution of about 2 At /t ~ 50%. A 
comparable data set for polyethylene on Vesuvio (Ref. 
[33h has been analyzed and the energy widths obtained 
are all substantially larger than the widths shown in Fig. 
[TBI except at the lowest energies where the resolution in 
the MARI experiment could easily be improved by us- 
ing a lower incident energy. When, however, the experi- 
mental widths are corrected for the Waller- Froman factor 
(Ref. H3) this decreases the width in the indirect geom- 
etry experiment by a large factor especially at the high- 
est scattering angles. The minimum energy widths of 
both experiments are then very similar. Therefore unlike 
previous results, our data cannot be reconciled with the 
indirect spectrometer data based upon the experiments 
having a different energy resolution. 

For Q=100 A -1 , experiments on polyethylene on Vesu- 
vio have obtained an energy width of ~ 25 eV where we 
measure ~ 7 eV. If we divide by the Waller-Froman Ja- 
cobain discussed earlier, the VESUVIO hydrogen width 
becomes 5 eV, comparable with the MARI results pre- 
sented here. Therefore, depending on how the data is 
analysed, the hydrogen recoil widths have been measured 
here to be either better or at least comparable to stud- 
ies on Vesuvio. We note that corrections for non Gaus- 
sian line shapes were found not to alter the conclusions 
obtained fron indirect geometry spectrometers^ Unlike 
previous results, we conclude that our data cannot be 
reconciled with the indirect spectrometer data based on 
a differing energy resolution. 

Through the use of several different incident energies 
with very different resolutions, we have shown that we 
obtain consistent answers both for the absolute cross sec- 
tion for hydrogen and for the constant momentum depen- 
dence. Therefore, the results do not depend on energy 
resolution as claimed in Refs. [31,32. The experimental 
resolution would be important if there was some break- 
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down in the Born Oppenheimer approximation when the 
line shape would be distorted by interaction with the elec- 
tronic energy levels. We suggest that these are unlikely to 
be observed and to cause strong deviations except where 
there is a strong coupling to a low energy mode possi- 
bly of order 100-500 meV or less. We have not observed 
any deviations in the line shape from polyethylene even 
though several electronic calculations have found energy 
levels with separation ~ 10 eV(Refs. 40.41) which are in 
good agreement with photoelectron measurements (Ref. 

We therefore conclude that the discrepancy in results 
from this experiment and the previous data based on in- 
direct geometry machines is due to experimental error. 
While this point needs to investigated through a series of 
careful experimental checks, Ref. [33j has suggested that 
understanding the flux distribution of the incident beam 
maybe one possibility for the problem. Another possi- 
bility could be the background introduced through the 
presence of a large amount of high-energy radiation from 
the initial interaction of the proton pulse with the tung- 
sten target. In conventional experiments on direct geom- 
etry instruments such as MARI, MAPS, and MERLIN, 
these high-energy particles are removed through the use 
of a nimonic chopper. Such a chopper does not exist on 
indirect geometry instruments and was not used in this 
experiment to access the high energy spectrum. While 
the background at short times is surprisingly large, we 
have found an excellent way to subtract this background 
that is comparable to the signal at incident energies of at 
least up to 100 eV. Presumably this background is also 
present on the indirect spectrometers but it is impossible 
to subtract this background on these machines. Such a 
background may not have been as severe in the case of 
the electron linac experiments discussed in Ref. [29ll47l . 
This point will need to be investigated in future work. 

While our results show conclusively the lack of a deficit 
and anomalous cross section, the experiments do show 
that direct geometry instruments can be used for large 
energy transfers. Previously, the use of direct geometry 
have been restricted to energies below ~ 500 meV (see 
for example Refs. I43ll44ll45l ) to study magnetic and elec- 
tronic excitations. There has been a considerable amount 
of interest in excitations at energy transfer greater than 
~ 1 eV. These have mostly been pursued with inelastic 
x-rays and light scattering where the condition of large 
energy transfers and small momentum transfers can be 
simultaneously achieved . 46 i 47 i 48 It will be a topic of fu- 
ture study to investigate whether neutrons can be used 
to study high energy magnetic and electronic excitations 
at energy transfers greater than ~ 1 eV. 

In conclusion, we have used a direct geometry chopper 
spectrometer to measure the recoil cross section of hy- 
drogen in polyethylene. We find the cross section to be 
independent of momentum transfer, in perfect agreement 
with values measured by Sears, and the results to be in- 
dependent of energy resolution. The data are in good 
agreement with theory based on the impulse approxima- 



tion. All of these results are in contrast to previous ex- 
periments on indirect geometry spectrometers and with 
theories based on quantum entanglement. We suggest 
that the discrepancy is the result of experimental issues 
using indirect geometry spectrometers. Our results do 
show the use and feasibility of direct geometry spectrom- 
eters at very large energies and we suggest future possible 
uses for these spectrometers. 
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VII. APPENDIX A: ENERGY CALIBRATION 

For energies where MARI is typically used (below ~ 
500 meV), the energy is calibrated from the time the 
neutron pulse reaches Mon 2 located just after the Fermi 
chopper (Fig. Because the monitors are sensitive to 
7 radiation (which are normally removed by the nimonic 
chopper) the short time spectrum of the monitors could 
not be used and therefore this method could not be relied 
upon for calibrating high incident energies beyond ~ 40 
eV. We therefore relied upon Mon 3 which is located 
further down the beam path implying it views less of 
the initial 7 flash from the protons striking the target 
and also the time of arrival of the neutrons will occur 
at considerably later times than the initial high energy 
background. Therefore, for our energy calibration and 
the low-angle detectors to calibrate the incident energy. 

Because the detectors and the monitors are based on 
different electronics, it is important to calibrate them 
against one another to ensure there are no systematic 
time differences. Fig. [151 a) illustrates a plot of the time 
of arrival of the neutron pulse at Monl as a function of 
time of arrival at Mon3. The data only includes spec- 
trum from incident energies less than ~ 40 eV as the 
time-of-flight peak position in Mon2 becomes difficult to 
reliably determine due to the background caused by the 
7 radiation normally removed by the nimonic chopper. 
The slope is the measured path length difference between 
the two monitors (Ljvfon3/ 1 Lmou2) and describes the data 
very well implying that Mon3 can be used to calibrate 
the energy of the incoming neutron pulse for higher inci- 
dent energies. This curve illustrates that the calibration 
on Mon3 is consistent and can be reliably compared with 
Mon2 at high incident energies or short times of flight. 

To test the electronics of the detectors and to ensure 
the time-of-flight of arrival is consistent for both He 3 de- 
tectors and glass bead monitors, Fig. [TS] b) displays the 
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FIG. 18: a) The time of arrival of the neutron pulse at Mon2 
as a function of time of arrival at Mon3. b) The time of arrival 
of the neutron at the low-angle detectors spanning 2#=[5°,7°] 
as a function of time of arrival at Mon3. 



time of arrival at the low-angle detectors with a vana- 
dium sample integrating over detectors ranging spanning 
28= [5° ,7°] . The detectors are located at sufficiently small 
momentum transfer that we expect that (for the case 
of vanadium) that no measurable recoil energy will be 
observable and the scattering will be dominated by the 
elastic incoherent cross section of vanadium. The fit in 
Fig. [15] b) is to a line with a slope equal to the ratio of 
the path lengths (Ldetl i^Mora) and with an intercept of 
3.5 \is. The intercept implies that the detector electron- 
ics count later than the monitors and the time difference 
must be corrected in calibrating the energy of the neu- 
trons, particular at high energies where the experiment 
occurs at very short times (~ 100-200 us). Based on the 
data presented in Fig. [THl we are able to calibrate the 
incident neutron energy at high-energies (or very short 
times) and also the energy transfer of the detected neu- 
trons. We hope that these studies will be useful in any 
future high-energy experiments on direct geometry spec- 
trometers. 
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FIG. 19: a) The full-width in energy squared (2r 2 ) of the 
hydrogen recoil line is plotted as a function of the final energy 
cubed (E 3 f ) for Ei=100 eV (panel a) and E;=40 eV (panel 6). 
The fit and parameters are described in the text. 



VIII. APPENDIX B: EXPERIMENTAL 
RESOLUTION 

The line-widths of the hydrogen recoil lines in energy 
are shown in Fig. [T2] for the different incident neutron 
energies and how they vary over the range of angles which 
we expect to be valid in our experiment. An unexpected 
feature of these results is the marked decrease in the line- 
width with scattering angle for large incident energies (a 
factor of ~ 3.5 for an Ei=100 eV). This result is sur- 
prising because, as discussed above, the width of the hy- 
drogen scattering is expected to increase with increasing 
angle in line with the impulse approximation whereas 
the experimental results show a steady decrease for large 
incident energies and large momentum transfers. Fur- 
thermore, it is commonly expected that the neutrons are 
produced in a very short time, about 0.4 /isec and the 
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chopper then takes a slice of incident energies by being 
open for a considerably longer time. The uncertainty in 
the incident energy (determined by the Fermi chopper 
opening) is then expected to dominate the energy res- 
olution for all energies. This is clearly not the case in 
this experiment especially for incident neutron energies 
of 100 eV and we will now consider this contradiction in 
more detail. 

In principle there are two substantial contributions to 
the energy width of a direct geometry time-of-flight spec- 
trometer. The neutrons are emitted by the moderator in 
a short pulse hrst which is about 0.5 fisec long and this 
will be treated as instantaneous. If the chopper is open 
for a time t then it will allow through all those neutrons 
which have the appropriate speeds to be able to pass 
through the chopper. These neutrons will have a range 
of velocities and that is given in terms of the energy by, 
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where E Q is the incident neutron energy, and L Q is the 
distance between the moderator and the chopper. This 
becomes 



AE 

E 



= 0.0027to 



(11) 



when the incident energy is expressed in eV, the time 
t in fisec and the distance in meters for the Mari di- 
rect geometry spectrometer. This energy width is then 
almost independent of the scattering angle apart from 
some smaller factors such as the thickness of the sample. 

The other effect that arises is that due to the neutrons 
having different velocities because the slower neutrons 
will arrive at the sample position after the faster neu- 
trons. This implies that in addition to the width coming 
from the different energies there is an additional width 
due to the spreading of the pulse with time. We shall as- 
sume that the spreading in time is given by f. A simple 
argument to estimate this contribution to the width is to 
estimate it as the opening time of the chopper multiplied 
by the distance from the moderator to the sample divided 
by the distance of the chopper from the moderator. This 
gives f = 1.17r . The neutrons have a definite time delay 
at the target which they continue to have throughout the 
instrument. However, the neutrons are recorded in time 
channels and to compare with the predictions in Fig. [T3] 
the spectra must be converted into energy spectra. We 
therefore relate the lengths of the neutron pulse in time 
to the final energy Ei by the following transformation: 
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2f 2Ei 



(12) 



When the energies are in eV, the length in metres and 
the time in /isec this becomes 
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(13) 



Assuming that these two contributions to the width of 
the scattering are independent we obtain the total width 
(AE) at the detector as 



AE= JAEl + AEl 



(14) 



This analysis can be tested against the experiment by 
plotting the full-width (2r) against the final neutron en- 
ergy Ei. In particular, a plot of (2r) 2 as a function of 
Ef gives the following expression for a straight line 



(2r) 2 = a + (3Ef. 



(15) 



with a = (0.0027r o ) 2 £;3 an d j3 = (0.0067f) 2 . The inter- 
cept and slope therefore provide a measure of r and f 
respectively. The data and fit plotted in this manner are 
presented in Fig. [T9l 

This expression gives, at least approximately, an ex- 
cellent description of the experimental data with 100 eV 
incident neutrons if f — 10 fisec and r =6 fisec. It also 
correctly describes the data taken with Ej=40 eV with 
t= 8 fisec and t d =11 fisec. The times fitted from the 
100 eV data have a ratio of 1.67 which is somewhat larger 
than predicted by the simple model described above. The 
data at 40 eV have a ratio which is smaller and the dis- 
crepancy indicates there are other contributing factors 
to the broadening of the time pulse. Possibilities include 
sample size, increasing transmission of the boron blades 
of the chopper with increasing incident energy, and error 
in the phase of the Fermi chopper. Despite these factors, 
the angular dependent widths are well described by our 
simple model governed by two time scales. 

We consider that this qualitative explanation of the 
resolution function is correct and that a full calculation 
including the size and shape of the scattering elements 
and the thickness of the detectors as well as dropping 
the independence of the time and energy aspects of the 
resolution would produce better answers. However, we 
do not wish to follow this here as we are confident that 
we have identified the main components of the resolution. 
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